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Abstract 

We discuss the possibility of the cosmic coiiiciciciicc generating the ratio of baryon asymmetry 
to dark matter in a Stueckelberg U{1) extension of the standard model and of the minimal su- 
persymmetric standard model. For the f7(l) we choose — Lr which is anomaly free and can be 
gauged. The dark matter candidate arising from this extension is a singlet of the standard model 
gauge group but is charged under L^ — Lt. Solutions to the Boltzmann equations for relics in the 
presence of asymmetric dark matter are discussed. It is shown that the ratio of the baryon asym- 
metry to dark matter consistent with the current WMAP data, i.e., the cosmic coincidence, can 
be successfully explained in this model with the depletion of the symmetric component of dark 
matter from resonant annihilation via the Stueckelberg gauge boson. For the extended MSSM 
model it is shown that one has a two component dark matter picture with asymmetric dark mat- 
ter being the dominant component and the neutralino being the subdominant component (i.e., 
with relic density a small fraction of the WMAP cold dark matter value) . Remarkably, the sub- 
dominant component can be detected in direct detection experiments such as SuperCDMS and 
XENON- 100. Further, it is shown that the class of Stueckelberg models with a gauged — 
will produce a dramatic signature at a muon collider with the ij(jLt+/i~ jX^ ,t^t~) showing 
a detectable Z' resonance while a{ii^ ji" e+e~) is devoid of this resonance. Within the above 
frameworks we discuss several broad classes of models both above and below the electroweak 
phase transition temperature. Asymmetric dark matter arising from a U{1)b-l Stueckelberg 
extension is also briefly discussed. Finally, in the models we propose the asymmetric dark matter 
does not oscillate and there is no danger of it being washed out from oscillations. 
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1 Introduction 



One of the outstanding puzzles in particle physics and cosmology relates to the so called cosmic 
coincidence, i.e., the apparent closeness of the amount of baryon asymmetry to the amount of dark 
matter in the Universe. Thus the WMAP-7 result, with RECFAST version 1.5 to calculate the 
recombination history jl], gives the baryonic relic density to be WOQEh^ = 2.255 it 0.054 and the 
dark matter relic density to be r^DM^o ~ 0.1126 ± 0.0036, which leads to 



4.99 ±0.20. 1 



The closeness of f^DM^o ^^^'^ ^B^o points to the possibility that the baryonic matter and dark matter 
may have a common origin; a possibility that has been noted for some time [2]. In this manuscript we 
analyze this issue in the framework of a Stueckelberg U{1) extension of the standard model (SM) as 
well as a Stueckelberg U{1) extension of the minimal supersymmetric standard model (MSSM) [3|-[9]. 
There are two main constraints in building models with asymmetric dark matter (AsyDM). First, 
we need a mechanism for transferring a B — L asymmetry produced in the early universe to dark 
matter. Second, we must have a mechanism for depleting the symmetric component of dark matter 
generated via thermal processes. 



The above issues have been discussed in the literature in a variety of works 10-19] (for a review 
see |20]). The models based on the Stueckelberg extensions we discuss in this work are different 
from the ones considered previously both in terms of the mechanism for depletion of the symmetric 
component of dark matter as well as regarding the implications for dark matter and signatures at 
colliders. Specifically, we consider a U{l)x extension of the standard model gauge group which is 
anomaly free. Further, we consider dark matter candidates which will carry lepton number but not a 
baryon number, and are singlets of the standard model gauge group. Now in the leptonic sector it is 
known [2l] that for the standard model case we may choose one of the linear combinations Lg — L^, 
Lfj, — Lr, Le — Lr to be anomaly free and can be gauged. The gauged Lg — has been discussed 
previously in the context of PAMELA positron excess and multi-component dark matter jsj and 



— Lj- in the context of muon anomalous moment 22 and in the context the PAMELA positron 
excess (23]. Here we consider a gauged — L-j- in the discussion of asymmetric dark matter as 
this choice is the more appropriate one for the analysis here. Specifically, we will consider a U{l)x, 
X = Lf^ — Lr Stueckelberg extension of the standard model as well as of the MSSM. As is well- 
known, the MSSM supplemented by supergravity soft breaking gives the neutralino as the lowest 
supersymmetric particle and with R parity a candidate for dark matter. Thus for the AsyDM to 
work in the MSSM extensions it is necessary to have the neutralino as a subdominant component. 
This issue will be addressed as well as the question if such a subdominant component may still be 
detectable in experiments for the direct detection of dark matter. It is found that the Stueckel- 
berg models with a gauged — Lr can produce dramatic signatures at a muon collider leading to 
a detectable Z' resonance in the fJ-~^n~ — )• fi'^n~,T~^T~ cross section while the H~^IJ.~ — >■ e^e~ cross 
section exhibits no such resonance. Finally we consider the possibility of an asymmetric dark matter 
in the Stueckelberg extension oi B — L. 

An important issue regarding asymmetric dark matter concerns the possibility that such matter 
can undergo oscillations [19]. Thus, for example, consider a model which allows for a Majorana 
mass term C = —niMXX + /i.e., where X is the dark particle. The presence of such a term along 
with other mass terms allows for the oscillation of X to its anti-particle X. Detailed analysis show 
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that in this circumstance over the age of the Universe the asymmetric dark matter would produce 
a symmetric component which would lead to pair annihilation. Such processes could completely 
wipe out the asymmetric dark matter generated in the early universe, and could render such models 



largely invalid 19 . In the models we consider mass terms that can generate oscillations are for- 
bidden because of gauge invariance, either U{1)l^-Lt or U{1)b~l- Thus dark matter oscillations 
are absent in the class of models we consider. However, we wish to add a further explanation here. 
In general there are two alternative possibilities for the Stueckelberg mechanism to arise. One pos- 
sibility is that it is a low energy remnant of a Higgs mechanism which contains a scalar field S 
that would in general allow a term of the type XX S^^ and thus generate a Majorana mass term 
for X as a consequence of VEV formation of the scalar field S. This would lead to dark matter 
oscillations once again. The second possibility, which is the view point adopted in this work, is that 
the Stueckelberg mechanism arises from the Green-Schwarz term in string theory (see, e.g., j24]). In 
this case there is no fundamental Higgs field which develops a VEV and thus Majorana mass term 
would not be generated and there would be no dark matter oscillation. 



We will discuss six broad classes of models. Three of these will be anchored in extensions of 
the standard model, one in extension of the two Higgs doublet model, and two in extensions of 
the minimal supersymmetric standard model. We will consider cases where the asymmetry transfer 
interaction may lie above or below the electroweak phase transition scale, i.e., the scale where the 
Higgs boson gets its VEV. In the supersymmetric case we will consider the case where all of the 
sparticles are in the thermal bath at temperatures where the asymmetry transfer takes place as well 
as the case where the first two generations of squarks are heavy and are Boltzmann suppressed in 
the thermal bath. 



The outline of the rest of the paper is as follows: In Section [2] we give a brief introduction to 
cosmic coincidence as well as asymmetric dark matter and the technique for the computation of the 
ratio of dark matter density to the baryonic matter density in the Universe and of the dark matter 
mass. Here we describe six broad classes of models which we will discuss in detail later. In Section [3] 
we carry out an explicit computation of these quantities in extensions of the standard model and of 
the two Higgs doublet model. In this section we also consider the case with inclusion of right-handed 
neutrinos. In Section |4] the analysis is redone for models in extensions of MSSM. In Appendix [A] we 
give a master formula for the computation of the asymmetric dark matter mass which is valid for 
temperatures above the electroweak phase transition scale. Here we show that the results of models 
discussed in previous sections can be deduced as limiting cases. In Section [5] we consider an explicit 
Stueckelberg extension of the standard model which generates asymmetric dark matter. We discuss 
solutions to the Boltzmann equations for relics in the presence of asymmetric dark matter, and show 
that the symmetric component of dark matter can be depleted from resonant annihilation via the Z' 
pole. In Section |6] we give a Stueckelberg extension of MSSM. Here a similar resonant annihilation 
of the symmetric component of dark matter is valid. There are several additional particles that 
arise in this case which include an extra scalar particle (the p) from the Stueckelberg sector. The 
decay width of this particle is computed in Appendix [B] and it is shown that it decays rapidly and 
is removed from the relativistic plasma. There are also additional neutralinos which we assume 
lie above the lightest MSSM neutralino and thus the lightest MSSM neutralino continues to be 
the lightest supersymmetric particle (LSP). In Section [7] we show that the MSSM neutralino is a 
subdominant component and thus does not interfere with the AsyDM mechanism. It is also shown 
here that the subdominant component can produce a spin-independent cross section which lies 
within reach of future experiments for the direct detection of dark matter. In Section |8] we discuss 
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the signatures of the models at a muon coUider. A [B — L) Stueckelberg extension is discussed in 
Section [9] which also produces asymmetric dark matter. Conclusions are given in Section 10 In 



Appendix [Cj we compute the /i''~/U — )• e^e at the loop level via Z' — 7, Z' — Z exchange and show 
that the corresponding production cross section is too small to be discernible. 

2 Cosmic coincidence and asymmetric dark matter 

In the analysis here we assume that a B — L asymmetry has been generated in the early universe. 
We do not speculate on how this asymmetry comes about as it could be by any number of different 
processes such as baryogenesis or leptogenesis [25| . Thus, for example, baryon asymmetry, specif- 
ically a non-vanishing B — L, can arise in the early universe by decay of super-heavy particles in 
some grand unified models |26[[27) consistent with experimental proton decay limits [28 1. A {B — L) 
asymmetry of this type will not be washed out by sphaleron processes which preserve B — L. Such 
an asymmetry is then transferred to the dark sector at high temperatures via an interaction of the 



form 12 



-^asy — -T-j—OriuOasy , (2) 

asy 

where Masy is the scale of this interaction^ and Oasy is an operator constructed from SM/MSSM 
fields which carries a non-vanishing B — L quantum number while Odm carries the opposite B — L 
quantum number. This interaction would decouple at some temperature greater than the dark 
matter mass. As the Universe cools, the dark matter asymmetry freezes on order of the baryon 
asymmetry, which explains the observed relation between baryon and dark matter densities. 

At the temperature where Eq. ([3| is operational, and using the fact that the chemical potential 
of particles and anti-particles are different, the asymmetry in the particle and antiparticle number 
densities is given by 
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^ gjT^ ^ {PlJ-iCi{h) bosons, 
6 \P^iiCi{f) fermions, 

where rii and fij denote the equilibrium number density of particle and antiparticle respectively, gi 



counts the degrees of freedom of the particle, Ei{q) = \J + mq where rrii is the mass of particle i, 
|Uj is the chemical potential of the particle (— jUj is the chemical potential of the antiparticle), and 
-|-1 (— 1) in the denominator is for the case when the particle is a fermion (boson). In the ultra 
relativistic limit (T ^ rrii) the mass of the particle can be dropped. In our analysis below we use 
the approximation of a weakly interacting plasma where /3/ii ^ 1, and /3 = 1/T and one has 

^ __ _ ^ hPfii + 0{{(3f,iy^) bosons, 

^ 6 ^ 1 /3^i + 0((/3^,)^) fermions. 



^ In the radiation-dominated era, the Hubble expansion rate is given by ~ /Mpi, where Afpi — 2.435x 10^* GeV 
is the reduced Plank mass. For an interaction suppressed by a factor 1/M"sy, the interaction rate at temperature T 
is r(r) ~ r^"+VA^asy- Thus, the interaction will decouple HT < H, i.e., when 



'asy 

-"asy 



Mf^y > MpiT^"-i . (3) 
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In the limit where Eq. (|5| holds we have Ci{b) = 2, Cj(/) = 1. This limit is a useful approximation as 
it simplifies the analysis of the chemical potentials that are needed in the generation of dark matter. 
However, full analysis can be easily done by using the exact expression of Eq. Q. We will discuss 
inclusion of these in Appendix [A| There the functions Cj(6) and Ci(/) introduced in Eq. Q will be 
found useful. The mass of the dark matter is constrained by the experimental ratio of dark matter 
to baryonic matter given in Eq. ([T]). Defining B to be the total baryon number in the Universe and 
X to be the total dark matter number, we obtain 

- 5, (6) 



f^mattcr B ■ UlB 

SO that the dark particle mass is given by 



moM ~ 5 • ^ • 1 GeV . (7) 



Applying the general thermal equilibrium method |29| (see also |11|), it is not difficult to express B 
and X in terms of the chemical potentials and then find their ratio. We note one subtlety is that 
while X and B — L (where by i? — L we mean the B — L in the standard model sector) are conserved 
after the interaction in Eq. ^ decouples, B is not. Thus, for example, the top quark would drop 
out from the thermal bath at some temperature Tf and one must solve the new set of ^ equations 
at T < Tt which would affect the computation of B although B — L is conserved. Typically one 
takes Tt to be Mf but it could be somewhat lower. Specifically, as the temperature drops below 
Mt ~ 173 GeV, the top quark becomes semi-relativistic but could still be involved in the thermal 
equilibrium constraints. A precise determination of Tj is out of the scope of this paper, and here 
we simply assume that Tt lies below Mt- Further, as the temperature falls below the temperature 
where sphaleron processes decouple, B and L would be separately conserved down to the current 
temperatures. Thus the relevant B to compute the dark matter mass in Eq. ([7| would be the baryon 
number below the sphaleron temperature which we label -Bfinal- It is useful to express X and -Bfinal 
in terms of B — L so that X = x{B — L) and -Bgnal = b{B — L) where b is to be determined later 
(see Eqs. (28) and (61)). Thus, Eq. ([T]) can be rewritten as 



moM ~ 5 • ^ • 1 GeV . (8) 

X 



We will discuss six broad classes of models labeled Models A-F (see Table[T]). Models A,B,C are 
anchored in the standard model while Model D is a two Higgs doublet (2HD) model. For Models 
A and D, the asymmetry transfer interaction, of the form of Eq. (|2]), is active only above the elec- 
troweak phase transition (EWPT) scale, i.e.. Tint > ^EWPT (^EWPT ~ 200 — 300 GeV where the 
Higgs gets its VEV). For Model B and C, the interaction which transfers the asymmetry could be 
active also below the EWPT scale, i.e., Tewpt > Tint- More specifically, in Model B we consider 
the temperature regime Tewpt > T^nt > Mt, and in Model C we discuss Tt > Tint > M^r {M\\r is 
the mass of W boson). Similarly, we discuss models based on extensions of the MSSM. Here we will 
focus on two cases; one of which is when Tint > -^SUSY (Model E) where MsusY is the (largest) 
soft breaking mass. In this case all the sparticles will be in the plasma. The second case (Model F) 
corresponds to when the first two generations of sparticles (with mass Mi) are heavy and drop out 
of the plasma (at some temperature Ti < Mi) while the third generation sparticles, the gauginos, 
the Higgses and the Higgsinos (with mass M2 <^ Mi) remain in the plasma. Thus for this case we 
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Model A 


SM 


Tint > TewPT 


Model B 


Tewpt > Tint > Mt 


Model C 


Tt > Tint > Mw 


Model D 


2HD 


Tint > TewPT 


Model E 


MSSM 


Tint > MsuSY 


Model F 


Ti > Tint > M2 > Tewpt 



Table 1: A list of six models which allow for generation of asymmetric dark matter. Models A,B,C 
are within the framework of the extensions of the standard model (SM) while Model D is an extension 
of a two Higgs doublet model (2HD). Models E and F are in the framework of an extension of the 
minimal supersymmetric standard model (MSSM). 



have Ti > Tint > M2 > Tewpt- 

These six cases are summarized in Table [T] There can be additional subcases for these models 
corresponding to different choices of the B — L transfer in Eq. 

3 Analysis in non-supersymmetric framework 

In this section we will determine the dark matter mass in terms of the B — L asymmetry in the 
non-supersymmetric framework utilizing Eq. ([s]) for Models A-D. We will discuss three different 
temperature regimes where the B — L transfer takes place and then deduce a general formula for 
computing the asymmetric dark matter mass. We note that the dark matter mass depends only on 
the {B — L)-charge of the operator Oasy that enters in Eq. ([2| and not on other particulars of the 
interaction. We will give several examples of the operator Oasy and compute the dark matter mass 
for them. 

3.1 T > Tewpt 

First we consider the case when the temperature is above the electroweak phase transition scale. In 
this case the following fields are in the relativistic plasma in the early universe: three generations 
of left-handed lepton doublets Li and quark doublets qi , three generations of right-handed charged 
leptons Cj and up and down- type quarks Uj and di {i = 1,2,3), and number \h of complex Higgs 
doublets Hi = {hf ,h^)'^. Since the Z boson and the photon couple to particle and anti-particle pairs 
they have a vanishing chemical potential. Further, in this temperature regime, SU{2)l symmetry is 
unbroken, the W and Z are part of the same gauge multiplet which requires that the chemical poten- 
tial of the W vanishes. The chemical potential of the gluon is zero and different color quarks carry 
the same chemical potential. The flavor (CKM) mixing among quarks ensures that the chemical 
potential of quarks in different generations are equal. However for the lepton sector, there is no such 



flavor mixing in the absence of neutrino masses 30 . Thus each of the lepton numbers {L^, L^^, 



for the three generations are separately conserved. Our notation is as follows: ^^.^^ei denote the 
chemical potentials of left-handed and right-handed leptons while /Xg- , ^u^ . , /i^ . stand for the chemical 
potential of left-handed and right-handed quarks. We assume that the chemical potential of all 
generations is the same and thus drop the subscript i and use for the chemical potential of the 
Higgs doublets (we assume all the Higgs doublets have identical chemical potential). 
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The Yukawa couplings 



^Yukawa = QeM^^^i + Ou^qiH^Ui + gd^iHdi (9) 

yield the following relations among the chemical potentials 

fJ-H = fJ-L - fJ-e = fJ-q - fJ-d = ^J'U - fJ-q ■ (10) 

Sphaleron processes {O^ph ~ nj=i 2 3 QiQiQi^i) S^^^ us one additional relation, 

3^g+^L=0. (11) 

The temperature at which sphaleron processes decouple is estimated to be Tsph ~ [80+54(m/i/120 GeV)] 
GeV |3l] . It is very likely that Tgph lies below TewpTj and thus the sphaleron processes are always 
active at T > Tewpt- Finally, the hypercharge neutrality condition requires the total hypercharge 
of the Universe to be zerc0 

^Hq + 6/Uu - - 3/UL - 3/Xe + '^^HfJ'H = . (12) 

Solving Eqs. ([9|-([l2) we can express all the chemical potentials in terms of the chemical potential 
of one single field, e.g., fj.^. Specifically one finds for Model A with Xh = 1 (suppressing a factor of 
/3TV6) 

5a = 3 X [2fig + {fin + ^^d)] = , (13) 

51 

La = 3 X {2nL + fJ-e) = y /UL , (14) 
so that {B — L)a = —-yI^l- And for Model D with \h = 2 we have 



15 , , 

Bd = -A^L , Lb = —yLL , (15) 



and {B - L)d = -f AiL- 
3.2 T < Tewpt 

Now we consider the case when the temperature is below the EWPT scale. After the Higgs gets its 
VEV, and the SU {2)i x U{1)y symmetry is broken, one has W^, Z, the photon and the Higgs scalar 
(/i) as the physical particles in the thermal bath. Again, since the Z and the photon only couple 
to two particles with opposite chemical potentials, their chemical potentials are zero. For tempera- 
tures above the top quark mass, the relativistic plasma includes three generations of left-handed and 
right-handed up-type and down-type quarks (uj^, um, dn and dj^), three generations of left-handed 



leptons {ciL and fj) and right-handed charged leptons {em), i = 1,2,3. As in Section 3.1 we will 
assume that the chemical potentials are generation independent. Thus dropping the generation in- 
dex we will use , , /^d^ , l^'d^ to denote the chemical potentials of left-handed and right-handed 
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The hypercharge of the Universe used in deducing Eq. ( 12 I is computed as follows 



y = 3 X [2 X 3 X l^i, + 3 X + 3 X + 2 x (-I)ml + (-2)/i.] + 2 x 2Ahmh , 

where the factor of 3 outside the first brace indicates summation over quark and lepton generations while inside the 
brace the factor of 3 for quarks indicates summing over colors, the factor of 2 for g, L and H counts two fields inside 
the doublets, and the additional factor of 2 for the Higgs is due to it being bosonic (see Eq. ([5|). 
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up- type and down-type quarks, /igj, and for left-handed leptons, for right-handed charged 
leptons, fiw for W~^, and fj^h for h. 

In the analysis below we make the following approximations: (1) At Tewpt > T > Mt, we still 
treat the top quark as relativistic gas. (2) AtTt>T> Mw, we treat the W boson as relativistic (all 
other particles, which have non- vanishing chemical potentials, are very light so the limit in Eq. ([s]) 
holds for them). (3) We assume Tt > T^ph, i-e., the top quark drops out of the thermal bath before 
the sphaleron processes decouple. 



For T < Tewpt ; the Yukawa couplings have the form 

>CYukawa = geALheiR + Qu^UiLhUiR + gd^dnhdiR + h.c. , (16) 

and since the Higgs boson is a real field and can couple to, for example, both eneiR and CiRCiL, we 
get 

= flh= fJ-UL - = f^dL - f^dR = fJ-SL - ^J■en ■ (17) 

Thus, the chemical potentials of left-handed and right-handed quarks/charged leptons are equal. 
The gauge interactions involving W bosons {C ~ W^fj^f) provide us the following relations, 

W = ^J'UL - l^dL ^UL+ di) , (18) 

IJ'W = fJ'u- IJ'eL {W'^ ^ + eix) . (19) 

The sphaleron processes give us one additional equation, 

^luL + +Hu = 0. (20) 

Since SU{2)l symmetry is broken below the EWPT scale, hypercharge is no longer a good quantum 
number. Further, the neutrality of the Universe now requires the total electrical charge to be zerc|^ 

2{fiuL + PtMfl + m) - {f^dL + f^dn + Ptei + /^e^) = . (22) 
Solving the new set of equations one finds for Model B 

36 

Bb=3x [{fi^^ +f^uR) + if^dL + fJ-dn)] = -y/^e , (23) 

75 

Lb = 3 X (/ie^ + fJ-en + fJ-u) = y Aie , (24) 
where we have expressed the results in terms of fi^ = /^e^ = ^e^j) and [B — L)b = —^fie- 

When the temperature drops below Tt, the top quark drops out from the thermal bath, and 
we are left with just five flavors of quarks. In this case (Tt > T > Mw) one must treat the first 



^ The result of Eq. ( 22 \ follows from the computation of the total charge Q which is given by 

(3 = 3 X [3 X + ^„^) + 3 X (-|)(Atd^ +P<ifl) + (-l)(Mei, +M.h)] + 2 X Sfiw, (21) 

where again, the factors of 3 for fermions outside the big brace indicates summing over generations, the other factor 
of 3 for quarks stands for summing over colors. For the W boson, 2 is the boson factor as given by Eq. ([sjl and 3 is 
the degrees of freedom of W. 
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two generations and the third generations separately. For the first two generations the analysis of 



Eqs. (17)-(20) still holds. For the remaining third generation leptons, we assume as before that 
the chemical potentials are identical to those for the first two generation leptons. Further, we note 
that the charge current process -f-)- ul + provides us with the relations nw = IJ-ul ~ M6l ^^^cl 
fif)j^ = ■ Thus we can treat Model C similar to Model B with only one modification to the charge 
neutrality condition, which now becomes 

4(/ini + /^«_r) + - 3{^J'dL + l^dn + fJ-CL + f^en) = . (25) 

Solving these equations we obtain for Model C 

90 

Be = 2{fiuL + ^J-ur) + 3(/id^ + ^dn) = -^l^e , (26) 



201 

Lc = 3 X (/iei + IJ-eR + f^u) = -J^fJ-e , (27) 



and {B—L)c = — ^/Ug. We note that the sphaleron processes will decouple below Tsph as mentioned 



already. Subsequently the baryon and lepton numbers would be separately conserved. Eqs. (17)- 



(19), and (25)-(27) would remain valid at Tsph > T > M\y- 



Following our assumptions given earlier, the top quark drops out of the thermal bath before 
sphaleron processes decouple. After the sphaleron processes decouple, B and L would be separately 
conserved. In other words, the ratio of B/{B — L) would freeze as soon as the sphaleron processes 
are no longer active. Thus, we obtain 

i^finai f B \ 30 _ _ ^28) 



B-L \B-L)n 97 
3.3 The AsyDM mass: non-SUSY case 

We discuss now in further detail the mechanism by which B — L \s transferred from the standard 
model sector to the dark matter sector and the determination of the dark matter mass. We con- 
sider the most general interaction which transfers the B — L asymmetry to dark matter at a high 
temperature: 

-CL^ = T^X'Ol^. , (29) 



asy jYjn asy 
-'"asy 



where the operator is constructed from the standard model fields, has a (B — L)-charge Q'^_i^, 
and X is the dark particle and has a (B — L)-charge Q^^^ = "Qij^'i/^l^ 

The parameterization of the asymmetric dark matter sector by the charge Q^^i is useful and 
we will utilize it in our analysis below. Also useful is the parameterization of the interactions in 
terms of the number of doublets and singlets that enter in Ofsy i i-^-i -^g> Nl,Nh numbers of q, L, H 
doublets and N^, N^, Ne numbers of uji, d^, singlets which are all active above the EWPT scale. 



Eq. (29) leads to the following constraints |13 



Nq^Jiq + Nl^IL + Nuflu + Ndfld + Nefle + NuflR + kfix = , (30) 

\Ng + \Nu + \Nd -NL-Ne + A;Qg^^ = , (31) 



The power of X can only be 2 or greater to ensure the stability of the asymmetric dark matter. 
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(32) 



Here Eq. ( 30 ) arises from the equilibrium of Eq. (29), Eq. (31 ) arises from the total (i3 — L)-cliarge 
conservation of the interaction, and Eq. (32) arises from the hypercharge conservation and the 
condition that the asymmetric dark matter must have zero hypercharge. Together with Eqs. (10)- 

11 



(IT2I), for Model A we obtain 
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If X is fermionic dark matter (EDM), we find, 

_ Xa 



79,, 



11 

'79 



Using Eqs. (pi) and (28), we obtain 



"T-FDM 



11.11 GeV 



Similarly, for Model D with two Higgs doublets, we have 



D _ 3 

-^Qb-lI^l 



so that 



W-FDM 



11.86 GeV 

^B-L 



(33) 
(34) 

(35) 

(36) 
(37) 



If the B — L transfer interaction is also active below the EWPT scale, the treatment is simi- 
lar. Assuming has N^, N^, N^, N,^, numbers of u,d,e,v,W~^ fields and recalling that at 
T < TewpT) the left-handed and right-handed quarks and charged leptons have the same chemical 
potentials, one finds the following constraints 



Nu^Xu + Nd^id + Ne^ie + N^fJ-u + Nw fJ-w + kfix = , (38) 



^N^ + \Nd-Ne-N, + kQ^^L = 0, (39) 



Nu-\Nd-Ne + Nw = 0. (40) 



We note that the last condition is from the charge neutrality of the operator O^^. Together with 
Eqs. ([17])- ([22]), we obtain for Model B, 

/^l = -yQg^L/^e. (41) 

The fermionic dark matter mass in this model reads 

B 15.60 GeV 

"^FDM ~ „0SM • v42j 

^B-L 

For Model C where the top quark is out of the thermal bath, we find 

^S = -§Qg^L/^e. (43) 
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and 

P 15.52 GeV 

"ipDM ~ ^^^SM — • (44) 

Now we consider the simplest example of the B — L transfer interaction {Q'^^j^ = —1) 

C^y = j^^l^^LH, (45) 

where "0 is the fermionic dark matter (which carries a lepton number of —1/3) and ^/^^ = ip'^'tp'ilj'^ . If 
this interaction is only active above the EWPT scale then the dark matter masses in Models A and 



D, and more appropriately in Models Ai and Di since the interaction of Eq. (45) is being used (see 



Table |2] which also includes a list of additional interactions) , are computed to be 

= 11.11 GeV Model Ai ; = 11.86 GeV Model Di . (46) 

If this interaction is also active below the EWPT scale, the dark matter masses in Models B and C 
are: 

= 15.60 GeV Model Bi ; = 15.52 GeV Model Ci . (47) 

Further, applying Eq. ([3| and the bounds in Table [T] one can estimate the mass scales for these 
interactions: 

> l-2x 10^ GeV, (48) 

1.2 X 10^ GeV > M^^ > 0.9 x 10^ GeV , (49) 

0.9 X 10^ GeV > Mfgi > 0.4 x 10^ GeV . (50) 



In the analysis above we focused on asymmetric fermionic dark matter. For bosonic dark matter, 

the masses would be half the fermionic ones, c.f., Eq. (p|. As an example, we consider now an 

iSM. 
asy ■ 



interaction with a higher dimensional operator Of^- 



C,sy = -^X\LHf. (51) 

-'"asy 

In this case, the dark matter could be either a fermion {X = ^, n = 4) or a boson {X = 4>, n = 3). 
This interaction gives rise to Models A2-D2 and Models A3-D3. As examples, for Models A2 and A3 



where Tint > ^ewpTj applying Eq. (35) we find that the dark matter masses are 



= 5.55 GeV Model A2 ; = 2.78 GeV Model A3 . (52) 



We explain now briefly the equality of asymmetric dark mass for the Models Ai,A4,A5,A6. 
From Eq. ^ we can write 

MOdm + /^oSNi = . (53) 
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1 Yk/nSM 



Model 



DM Mass 



Model 



DM Mass 



Model 



DM Mass Model 



DM Mass 



Ai 
A2 
A3 
A4 
As 
Ae 



11.11 GeV 
5.55 GeV 
2.78 GeV 
11.11 GeV 
11.11 GeV 
11.11 GeV 



Bi 
B2 
B3 
B4 
B5 
Be 



15.60 GeV 
7.80 GeV 
3.90 GeV 
15.60 GeV 
15.60 GeV 
15.60 GeV 



Ci 
C2 
C3 
C4 
C5 



15.52 GeV 
7.76 GeV 
3.88 GeV 
15.52 GeV 
15.52 GeV 
15.52 GeV 



Table 2: A display of the various interactions that allow a transfer of the B 
the standard model sector to the dark matter sector. 



Di 
D2 
D3 
D4 
D5 



11.86 GeV 
5.93 GeV 
2.96 GeV 
11.86 GeV 
11.86 GeV 
11.86 GeV 



L asymmetry from 



For Models Ai,A4-A6 we have 

LH (Ai) 
LLe^ (A4) 
Lqd'' (A5) 
u^d^d" (As) 



asy,l 

/i/nSM = 2fil — fie , 

asy,4 

flQSM ^ = fJ-L + fJ-q — fJ-d , 



asy , 6 



From the /i equations Eqs. (10) and (11), it is easy to see that 



asy,l asy, 4 asy, 5 



SM 

asy, 6 



(54) 
(55) 
(56) 
(57) 

(58) 



Eq. (58) implies that the dark matter has the same mass for the Models Ai,A4-A6. Similar analysis 
holds for Models Bi,B4-B6, Ci,C4-C6 and Di,D4-D6. 



We summarize all our results in Table [2] where we list the dark matter mass for the various 
interaction^ that can transfer the B — L asymmetry from the standard model sector to the dark 
matter sector. We note that for the first five interactions, the dark matter carries lepton number, 
while for the last one, it carries a baryon number. 



3.4 The AsyDM mass: including the right-handed neutrinos 

In the analysis above we used the framework of the standard model where we have no right-handed 
neutrinos and the neutrinos are assumed massless. The nature of neutrino masses is currently not 
known, i.e., whether they are Majorana or Dirac, but in the context of a gauged Lf^ — Lr symmetry it 
is more natural for the neutrinos to have Dirac masses which implies that we introduce right-handed 
neutrinos, one for each generation. We discuss now the effect of this inclusion on the analysis, i.e., 
on the fi equations, on the B /[B — L) ratio and thus on the DM mass. 



Since the right-handed neutrino u^i has hypercharge and electrical charge, it does not affect 
! neutrality conditions (such as Eq. (12) or Eq. (22)). i'r is only involved in one interaction 
^ LiH^um before the electroweak phase transition (or C ~ i^nhviR after EWPT), which gives us 



^ In the first column of Table[2] L, H and q stand for SU(2)l doublets as discussed in T > Tewpt regime (Model A 
and D). When the temperature drops below EWPT scale (Model B and C), since SU(2)l symmetry is broken, these 
interactions should be rewritten in terms of the contents of the original doublets. We omit this step for simplicity. 
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1 Yk/nSM 



Model 



DM Mass 



Model 



DM Mass 



Model 



DM Mass Model 



DM Mass 



M: 



a; 

A' 
A' 
A' 
A' 



12.12 GeV 
6.06 GeV 
3.03 GeV 
12.12 GeV 
12.12 GeV 
12.12 GeV 



B'4 

B'5 
B' 



15.58 GeV 
7.79 GeV 
3.90 GeV 
15.58 GeV 
15.58 GeV 
15.58 GeV 



C'4 



15.52 GeV 
7.76 GeV 
3.88 GeV 
15.52 GeV 
15.52 GeV 
15.52 GeV 



D2 

D4 
D' 



12.70 GeV 
6.35 GeV 
3.18 GeV 
12.70 GeV 
12.70 GeV 
12.70 GeV 



Table 3: A display of the various interactions that allow a transfer of the B — L asymmetry from 
the standard model sector (including the right-handed Dirac neutrinos) to the dark matter sector. 



= l-i-Uiii ~ f^Li before EWPT (or fi^.^ = fi^.^^ after EWPT). Thus the only change would be the 
total lepton number since carries lepton number 1. By including the right-handed neutrinos, a 
reanalysis gives the following formulas 



A' 
m-FDM 



"^FDM 



12.12 GeV 
15.58 GeV 



"^FDM 



"IpDM 



12.70 GeV 
15.52 GeV 



(59) 
(60) 



where we use a prime to denote all the models with the right-handed neutrinos. With the inclusion 
of right-handed neutrinos in the thermal bath 6' is determined to be (c.f. Eq. ([8|) 



b' 



^final 

B- L 



B 



B-L 



5 

21 



0.24. 



(61) 



With above formulas, we compute the dark matter masses for the same B — L transfer interactions 
displayed in Table [2] and they are collected in Table [3] We note that inclusion of right-handed 
neutrinos generates less than a 10% effect at most and no effect for Model Cj^ 



4 Analysis in supersymmetric framework 

We now consider the analysis in a supersymmetric framework specifically within an extended MSSM. 
Since the supersymmetric case can have its own dark matter candidate, i.e., the neutralino, the relic 
abundance of the neutralino must be depleted. For this reason, we only consider the parameter 
space where relic density of the neutralino is much smaller than the WMAP value for cold dark 
matter (CDM) and is thus only a subdominant component. Below we discuss two regimes, one 
where Tint > MsusY and the other where Ti > Tint > M2 > Tewpt- 

® The reason for this is simple, i.e., using Eq. Q and recalling the fact that Bfinai = Be = Be, since the inclusion 
of right-handed neutrinos does not change the total baryon number, we have 

mom = 5 • ■ 1 GeV = 5 • ^ • 1 GeV = 5 • ^ • 1 GeV = ■ 
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4.1 T> MsusY 

In this regime since the temperature is above the SUSY breaking scale all sparticle masses must be 
included in the equations. This case is very similar to the discussion of T > Tewpt in the standard 
model framework, except this time our particle spectrum includes all the standard model particles, 
the extra Higgses as well as the sparticles. For brevity we will use the same symbols for the chemical 
potentials, though now they stand for not only the standard model fields, but also their super- 
partners. The chemical potential equations obtained from Yukawa couplings and sphaleron processes 
remain the same. The only equation modified would be the hypercharge equation, which become^ 

3fig + 6;u„ - - - 3/Ue + 2fj.H = . (62) 

Solving the chemical potential equations, we find that for Model E the total baryon and lepton 
numbers are given by 

5e = 3 X 3 X [2ng + + fid)] = -12/iL , (63) 

153 

Le = 3 X 3 X {2fXL + /Ue) = — /iL , (64) 

so that (B — L)e = —^Y'fJ.L- Note that in the above equations, the extra factor of 3 = 1 + 2 (compare 
to the standard model case) takes into account the contributions of both fermions and bosons from 
the superfields, c.f. Eq. 

4.2 Ti > T > M2 > Tewpt 

Here we consider two soft breaking mass scales Mi and M2 where Mi ^ M2. When temperature 
drops below Ti, all the super-particles with masses greater than Mi would drop out of the thermal 
bath. We assume that this is the case for the first two generations of squarks and sleptons. Similar 
to Model C, we simply assume here that these super-particles would drop out of the thermal bath 
at Ml > Ti > M2. Thus the only super-particles remaining in the thermal bath are the third 
generation sparticles, the gauginos, the Higgses and the Higgsinos. We make the approximation 
that these particles are relativistic at Ti > T > M2. This case is labeled Model F. Following the 



analysis of Eq. (62) we find that the vanishing of the hypercharge for Model F gives 

+ 10/i^j - - 5/XL - bjjLe + 6/xj/ = . (65) 
Solving the ^-equations, we obtain 

20 

= (3 X 1 + 2) X [2^g + {ixu + ^^d)] = -y/^L , (66) 

485 

Lf = (3 X 1 + 2) X (2/iL + fie) = —fiL , (67) 

and (S-L)f = -If^L. 



^ The hypercharge of the Universe for the case when T > Tsusy is given by 

F = 3 X {3 X [2 X 3 X + 3 X + 3 X {~l)fid + 2 x (-I)ml + (-2)^=] + 2 x (fin^ - fi„^)}, 

where the counting is similar to discussion in footnote [2] The Higgs mixing term in the superpotential, i.e., W = 
nHuHd indicates + I^h^ = 0, and so we define hh = f^H^ = ^I^h^- 
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4.3 The AsyDM mass: SUSY case 



The supersymmetric interactions which transfer B — L asymmetry typically have a different form 
than the ones in the non-supersymmetric case. The most general interaction that transfers B — L 
to the dark sector for the MSSM case is 

W^asy = -^X^O^^^ , (68) 

where the dark matter superfield X = {(j)x,''Px) with (px as the bosonic and ipx as the fermionic 
component. Now the following possibilities arise in terms of dark matter. First, after soft breaking 
if (fix and ipx have a similar mass, both of them are stable, and could be dark matter candidates. 
Next, consider the case where one of the components has a much larger mass than the other and 
would decay into the lighter one. In this case we have two possibilities: either (px is heavier than 
ipx so that (px — ^ i^x + X^* (where is a Stueckelberg neutralino) in which case ipx is the dark 
matter candidate, or ipx is heavier than (px so that ipx — ^ (t>x + X^* in which case (px is the dark 
matter candidate (The possibility that either or the MSSM neutralino is a dark matter candidate 
is discussed in Section [6]). For either of these three cases, when computing the total dark particle 
number from Eq. (|5]), we need to multiply by an additional factor of 3, since both bosonic and 
fermionic components of the dark matter superfield would contribute. But for concreteness in our 
analysis we will assume that ipx is lighter than (px and thus would be the asymmetric dark matter. 



Applying the same method we used in Section 3.3 we find 



E 11.11 GeV F 6.51 GeV 



Thus for the B — L transfer interactions with Q^^^^"^ = —1, where O^y^'^ can be LH^, LLe^, Lqd'^, 
or u'^(F(F, the dark particle masses are 

mx = 11.11 GeV Model E ; mx = 6.51 GeV Model F . (70) 

the dark 



For the case Wasy = . A X'^{LHu)^ with = — 2, which we will discuss in Section 

particle masses are 



mx = 5.55 GeV Model E ; mx = 3.25 GeV Model F , (71) 
and using Eq. ([3| one finds 

>3.7x 10^ GeV. (72) 



If we include the right-handed neutrinos in the supersymmetric framework, the analysis is similar. 



and Eq. (69) becomes 



12.12 GeV 

^0MSSM 



F' 



6.99 GeV 

^B~L 



(73) 
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5 Asymmetric dark matter in a Stueckelberg extension of the SM 



As discussed in the Introduction, one of the major problems for an acceptable AsyDM model is 
to have an efficient mechanism for the annihilation of dark matter that is produced thermally. In 
general one has 

where is the relic density of asymmetric dark matter (which carries a nonzero (B — L)-charge) 
and is the relic density of dark matter which is produced thermally. For the asymmetric dark 

matter to be the dominant component, one must significantly deplete the symmetric component of 
dark matter. Specifically we will use the criteria that f^^^ / ^^dm < O.lj^ Thus we investigate if the 
symmetric component of dark matter produced by thermal processes can be annihilated efficiently. 
We accomplish this via the exchange of a gauge field using the Stueckelberg formalism where the 
gauge field couples to — Lr- 



For illustration let us consider Model Ai, which is governed by the interaction Eq. (45) operating 
at Tint > Tewpt- The corresponding dark matter mass is 11.11 GeV. Further, we require the dark 
matter particles ip to have a non-vanishing fi or r lepton number. The total Lagrangian is given by 

^ = ^SM + ^u{i) + -^St , (75) 

where >C(/(i) is the kinetic energy for the gauge field for the — Lj- symmetry, and for Cst we 
assume the following form: 

Cst = -liMcC^ + d^af . (76) 

In the unitary gauge the massive vector boson field will be called Z' and its interaction with fermions 
in the theory is given by 

Ant = IgcQti^i'^i'C^ + \gcQUi^fC^ , (77) 

where / runs over /x and r families and Q(j = —Qc- 

5.1 Resonant annihilation of symmetric dark matter 

We discuss now the details of the annihilation of the symmetric component of dark matter. We 
will show that the relic density for such dark matter can be reduced significantly below the WMAP 
value with resonant annihilation via the Z' pole, i.e., via the process ipij ^ Z' ^ //0 Thus, by 



using Eq. ( 77 ) one can compute the tpi/j — >■ // annihilation cross section and using the Breit-Wigner 



form for a resonance one has 



(T.,j,^,f = a^\{s- Ml, + iTz'Mz')\ , (78) 



* The analysis of previous sections was based on the assumption Sl^^ / Jlmattcr ~ 5. Inclusion of a small contribution 
(i.e., < 10%) of symmetric component to dark matter will proportionately affect the determination of the dark matter 
mass. It is straightforward to take account of this contribution but we do not carry it out explicitly as it is a relatively 
small effect. 

^ While the thermal dark matter can annihilate into second and third generation leptons at the tree-level, such 
an annihilation into the first generation leptons can come about only at the loop level involving the second and third 
generation leptonic loops. Thus the annihilation of thermal dark matter into first generation leptons is significantly 
suppressed relative to the annihilation into the second and the third generation leptons. 
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647rs/3^ 



Clip 



s\l + ^pjPl) + 4M2(s - 2m}) + Amjis + 2M^) 



where f3f^^ = (1 — Am} ^/s)^^'^ . The relevant partial Z' decay widths are given by 

/I .\2 

r{z' i^i)) = 



1 nfX „ 



:9cQ. 




AM':. 



2\ V2 



Ml 



G {Mz' - 2M, 



(79) 

(80) 
(81) 



where = 1 for f = fj,,T and rj = 1/2 for f = i^fiji^T- A constraint on (7c comes from the 
contribution of the Z' to — 2 32 33 , which is given by 



A(ff^ - 2) 



(82) 

In the analysis here we impose the constraint that the Z' boson contribution be less than the exper- 
imental (4o-) deviation of Aa^ = A((5^ - 2)/2) = (3.0 ± 0.8) x 10"^ |32|[33| , which is the constraint 
commonly adopted in analysis of supergravity based models. 

The relic densities of tp and are governed by the Boltzmann equations and have been discussed 



in several works 14 15 for the case of asymmetric dark matter. Typically it is seen that the effect 



of including the asymmetry in the Boltzmann equations lead to a significant effect on the relic 
density. In these works the analysis was done in the approximation {av) = a + bv'^. Our analysis 
below differs from these in that for our case annihilation via the Z' pole is the dominant process. 
Thus in our analysis we need to carry out an explicit thermal averaging over the Breit-Wigner 
pole. It is convenient to work with the Boltzmann equations for the quantities = n^/{hT^), and 



rJ{hT^) where 



■ ) is the number density of particle ip {ip) and the combination hT^ 



appears in the entropy per unit volume, i.e., s = (27r^/45)/iT^ where h is the entropy degrees of 
freedom. The Boltzmann equations obeyed by and take the form 



dx 
dx 



^ = «(-)(w^-/;vj'), 



where x = fc^T/m^ in which fc^ is the Boltzmann constant and hereafter we set ks 
given by 

^vr V Ag dry 



(83) 

(84) 
1, and a is 

(85) 



where g is the degrees of freedom that enter in the energy per unit volume, i.e., p = ^gT^, where 
T(t) = T^{t) is the photon temperature. Numerically a{T) = 6.7 x 10^*^ GeV^ for g = h = 68 at 
T = 0.5 GeV. {av) is the thermally averaged cross section 



{av) 



dv {av)v e 



2 -v^/ix 



dvv\-^"l^^ 



(86) 
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Further, in Eqs. (83) and (84) and are values of and at equilibrium. Now one can 
obtain the result from Eqs. (83) and (84) that the difference of and i.e., 



1 = U- h, 



(87) 



is a constant. Assuming that the asymmetric dark matter currently constitutes a fraction A of the 
dark matter relic density, one can evaluate 7 to be 



7 ~ A ■ 



5pc 



A70, 70 ~ 1.3 X 10"^" (m^ ~ 10 GeV) , 



where the 5/6 in 70 is due to Eq. ([T|. 

It is now straightforward to obtain the individual relic densities for ip and ip. Thus one integrates 
Eqs. (83) and (84) from the freeze-out temperature to the current temperature of X? = 2.73 K. In 



the integration we will make the following approximation which is conventionally done, i.e., we move 
a out of the integral and replace it with a(xj), i.e., by the value of a at the freeze-out temperature. 
The matter density of ip at current temperature is given by = m^n^{xo) where xo = T^/m^ and 
TO is the current photon temperature of 2.73 K. The relic density then is 



= m^n^{xo) / Pc 



(89) 



where pc is the critical matter density so that />c = (3 x 10 GeV) /iq where ho is the Hubble 



parameter. The integration of Eq. (83) straightforwardly gives 



2.2 X 10 



-11 



g{xf)h{xo,Xf) 



2.73 



where 



J{Xj 



(av) dx , h{xQ,Xf) 



hjxo) 
h{xf) 



fi,{xf) 



-iJ{xf) 



(90) 



1 + 



1 (Tdg 



4 \gdT 



(91) 



and ^ = 0(3:^)7 where a{xf) is the value of a evaluated at the freeze-out temperature, and where 
g{xf) {h{xf)) are the energy (entropy) degrees of freedom at freeze out and h{xo) is the entropy 
degrees of freedom at the current temperature. The derivative term \i'^^)xf is small and is often 
dropped, while h{xo) = 3.91 30,34 and we estimate h{xf) ~ g{xf) given Tj. As discussed below Xf 
is typically of size ~ 1/20 and thus Tj = m^Xf ~ 0.5 GeV for ~ 10 GeV. Now for Tf ~ 0.5 GeV, 
h{xf) ~ 68 which gives h{xo,Xf) ~ 1/17.5. The quantities f^{xf) {f^{xf)) are (/^) evaluated 
at freeze out. Analogous to the relic density for ip, we can get the relic density of ip by integration 
of ^ and we obtain 



n^h^ = 2.2 X 10 



-11 



g{xf)h{xo,Xf) 



The total dark matter relic density is 



2.73 



Ui^f) fJ(Xf) _ i 



(92) 



^DM = ^ 



(93) 



18 



From Eqs. (90) and (92) one obtains the ratio of the current rehc densities of and ip to be 



U{xf) 



(94) 



The front factor f^{x j)/ f^{xf) in Eq. (94) takes into account the asymmetry that exists at the 
freeze-out temperature. The size of this effect is estimated at the end of this section and could be 



as much as 20%, and thus significant. Our result of Eq. (94) is in agreement with the analysis of 14 



We discuss now the evaluation of the freeze-out temperature. The definition of this quantity 
differs somewhat in various works (see, for example, p5||36| ) but these differences are not very 
significant. We adopt here the definition of |35| where the freeze-out temperature Tf is defined as 
the temperature where the annihilation rate per unit volume equals the rate of change of the number 
density. For our case this implies 



dx 



^'\'^'^)f^f7 ' at x = X/ = Tf/m^ 



where /^"^ takes the form 



a^x-'l^e-'/^ 



(95) 



(96) 



where = g^{2Ti)-^l^h-\T) « 
freedom of the dark particle {g^ 
determined by the relation 



9.3 X 10 around T = 0.5 GeV, and denotes the degrees of 
= (7^ = 4 for Dirac spinors). The freeze-out temperature is then 



(97) 



L-j )^ ■ — a{av){a^ + ^x^pe^/^f) . 

For the case of no asymmetry, i.e., in the limit 7 — )• 0, Eq. (97) reduces down to the well-known 
result 35 . One may compare the analysis of the freeze-out temperature given by Eq. (97) with the 

(98) 



one using the alternate criterion |36 



A(x;) 



where A(x) 



(/^(x) — f^{x)) and c is order unity. Using Eq. (96) in Eq. (98) one gets 
(x}/^ - Ix'^'"^ - a7(at;)x3/2) e^l^i = aa^c{c + 2){av) . 



(99) 



For 7 = 0, Eq. (99) reduces to the result of [36| while 7 7^ gives the correction due to asymmetry. 
Further, we see that Eq. (99) reduces to Eq. (97) when c = — 1. To compute the sensitivity 
of the freeze-out temperature on the asymmetry it is useful to utilize the scale factor A defined in 
Eq. (88). On using Eq. (97) we can obtain an approximate expression for dxf/dX so that 



dx//dA 



(100) 



From above we can compute the first order correction to the freeze-out temperature due to the 
asymmetry. To the leading order one has 



1 



:S(xO)5/2ei/"/ 



:ioi) 
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where Xj is the zeroth order of the Xf, i.e., when 7 = 0. We note that the correction to the freeze-out 
temperature due to asymmetry is independent of (av) to leading order. Using = 3.7 x 10~^, 
Xf = 1/17.5 and 7 = 70 = 1-3 x 10"^*^, one finds that the correction to Xf is around a percent for 
the choice of the parameters given. Further, as 7 (and hence ^) — )• 0, one has — )• 1 and in this 
limit one has 

n^hl = n^hl = 2.2 X 10-'' ./^)h{xo,Xf)(^) (102) 



Now rapid annihilation of dark matter can occur if the sum of the dark matter masses is close to 



the Z' pole and there is a Breit-Wigner enhancement 37,38 . Thus for the case we are considering if 
the mass of the Z' is close to twice the mass of the dark particle, then one can get a large annihilation 
cross section and correspondingly a small relic density. An analysis of the relic density arising from 
the annihilation of symmetric dark matter is given in Fig. [T] and the analysis shows that the relic 
density arising from the symmetric component of dark matter can easily be made negligible, i.e., 
less than 10% of the cold dark matter density given by WMAP. In Fig. [T]we give the analysis for the 
case with no asymmetry, i.e., 7 = (left panel) and the case with asymmetry (right panel) where 
7 = 1.3 X 10^^'^. A comparison of the left and the right panels shows that inclusion of the asymmetry 
has a substantial effect on the relic density. Specifically it further helps deplete the relic density of 
ip (the symmetric component of dark matter). For the case of g'c' = 1 the allowed upper bound of 
the Z' mass increases by about ~ 100 GeV in the presence of an asymmetry when 7 = 1.3 x 10^^*^. 
It is also instructive to examine the ratio of the thermal relic density for the cases with and without 
asymmetry. Here one has 

_ (^^^g)7=70 ^ (Xf) , . 

As ^ — 7- 0, f'-{^J^^^ ~^ 1 s-iid thus i? — 7- 1. However, if we assume that the asymmetric dark matter 
is responsible for 5/6 of the total relic density, then for ~ 10 GeV, one has 7 = 1.3 x 10^^*^ 
and f{jj{xf) = 6.8 x 10~^^ which gives j.'^l^j^j = (1 + ^0/ f^{xf)) ~ 1-2. In this circumstance one 
finds that R is always less than 1. Thus one finds that the inclusion of asymmetry helps deplete the 
symmetric component of dark matter. 



6 Asymmetric dark matter in a Stueckelberg extension of the MSSM 

The analysis of dark matter in the MSSM extension is more complex in that there are now three 
contributions to the dark matter relic density, i.e., from the asymmetric and symmetric components 



as in Eq. (74) and from the neutralino. Thus here one has 

where O^o is the relic density from the neutralino. In this case for the asymmetric dark matter 
to work, one must significantly deplete not only the symmetric component of dark matter but also 
the contribution from the neutralino. Thus here we take the criterion that f^^^ / ^dm < 0.1, and 
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Figure 1: (color online) An exhibition of the thermal relic density of ^/^ as a function of Mz' in 
the model with gauged — for different values of the coupling constant. The left panel shows 
the case 7 = and the right panel shows 7 = 1.3 x 10"^''. In both cases, the analysis shows that 
the component of dark matter that is thermally produced can be efficiently depleted by resonant 
annihilation via the Z' pole. 



VL^a /r^DM < 0.1. For the analysis of AsyDM in extensions of MSSM we consider the interaction 

W^, = -^X\LH^f. (105) 

asy 

Here we note that the choice VFasy ~ X^LH^ would have allowed the decay — )• XXu ■ ■ ■ 
and would have required the constraint m^o < 2mx for the neutralino to be stable. Further, 
while the choices Wasy ~ X^LLe'^, X^LqdP do not allow the neutralino decay at the tree-level, 
such a decay can occur at the loop level since it is not forbidden by a symmetry. Additionally 
Wasy ~ X'^LHu, X'^LLe'^, X'^Lqd'^ can also preserve the stability of the neutralino. Here and else- 
where we are assuming that the Stueckelberg neutralinos are heavier than the lightest neutralino in 



the MSSM sector (see the discussion following Eq. (113)). Returning to Eq. (105), the corresponding 



dark particle masses are computed to be 5.55 GeV (Model E) and 3.25 GeV (Model F). Now the 
Stueckelberg extension of MSSM, is more complex than the SM extension. We exhibit the relevant 
parts of this extension below. 

For the Stueckelberg Lagrangian of the supersymmetric case we choose [S] 

Cst = J de^de^ [MC + s + s]\ (loe) 

where C is the U{l)c vector multiplet, S and S are chiral multiplets, and M is a mass parameter. 
We define C in the Wess-Zumino gauge as C = -Baf'dCf, + iOOeXc - imXc + leOMDc , while 
S= ^(p + ia)+ex + iOa^'md^,p + id^,a)+ OOF + ^OeOa^'df.x + lOeeeinp + ina) . Its complex scalar 
component contains the axionic pseudo-scalar a, which is the analogue of the real pseudo-scalar 
that appears in the non-supersymmetric version in jsj. We write Cst in component notation as (see 
e.g. D) 

-Cst = -^(MC^ + d^af - ^{d,,pf - ix<y^d^x + 2|i^|' + MpDc + Mx'Xc + MxXc • (107) 
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For the gauge fields we add the kinetic terms 



>Cgkin = -\c^,C^' - iXc^^d^Xc + \dI , (108) 

with C^jy = dfj_Cy — duC^. For the matter fields (quarks, leptons, Higgs scalars, plus hidden sector 
matter) chiral superfields with components {fi,Zi,Fi) are introduced and the matter Lagrangian is 
given by 

/:matt = -\D^Zi\^ - iha^Dji - (iV2gcQcZifi~Xc + h-c ) + gcDc{ziQcz^) + \Fi\^ , (109) 



where {Qc^Qc) ^re the charge operator and coupling constant of U{l)c, and = + igcQcC^i 
is the gauge covariant derivative. It is convenient to replace the two-component Weyl-spinors 
iX) X) ^ {Xc ,Xc) by four-component Majorana spinors, which we label as ips = (Xo)X")"^i ^-ncl 
Ac = (Aca, A^). The total Lagrangian of the MSSM then takes the form 



with 



-^^StMSSM = -^^MSSM + ^U{1) + ^-^St , (HO) 



-\i^sl^d^i^s - -^C^uO^" - '-Xcrd^Xc + MijsXc 



^Zi\^ - \D^Zi\'^ + pgcM{ziQcZi) 



i 

+\9cCjirQcfi + V2gc{iz^QcI^Xc + h.c. )] - \[9cY. ^^Qcz}^ ' . (Ill) 

i 

As in the SM case we assume that the C/(l)c is a gauged — L^-. Further, we assume that all 
hidden sector fields while charged under \J{\)c are neutral under the MSSM gauge group and some 
of the MSSM particles, i.e., the second and the third generation leptons, are charged under \J(\)c- 
As discussed already an essential ingredient to explain the cosmic coincidence is that the symmetric 
component of dark matter produced in thermal processes is significantly depleted. For the MSSM 
Stueckelberg extension the analysis of annihilation is essentially identical to the case of the Stueck- 
elberg extension of the standard model and we do not discuss it further. 

We now discuss the fate of the extra particles that arise in the f/(l)c Stueckelberg extension of 
MSSM. This extension involves the following set of particles: Z' , p,ip, (fyjips^ Xc- The decay of the 



Z' has already been discussed. Next we consider the p. Eq. (Ill ) gives the interaction of the p with 



the sfermions. Specifically its couplings to the mass diagonal sfermions are given by 



^pPI = -9pMp cos(20^.) ( flJu - fy2^ ) + sin(20^.) ( /jja. + flJu ) , (112) 



where ft refer to p,v^,T,Vr. Thus the p will decay via second and third generation slepton loops 
into f^z/^, r+r~, Vt^Ut and disappear in the thermal bath quickly (see Appendix |b]) . Next 

we discuss the neutralino sector. Here in the U{\)c Stueckelberg extension of MSSM the neutralino 
sector is enlarged in that one has two more fields, i.e., the gaugino, and the higgsino fields {^s^^c) 
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as mentioned earlier. In this case the neutrahno mass matrix of the U{l)c extension of MSSM is 
given by 



M 



ncutralino 



Mst 


02x4 


04x2 


-A^MSSM 



where Mst is in the basis (^'5',Ac), M is the Stueckelberg mass and M is the soft mass. The 
neutrahno mass eigenstates arising from Eq. ( 113 ) can be labeled Xi^, xf*- We consider the possibility 



that the Stueckelberg neutralinos are heavier than the LSP of the MSSM (x?) and decay into the 
MSSM neutrahno which is assumed to be stable. In this case one will have more than one dark 
matter particle, i.e., the ■0 from the Stueckelberg sector and Xi from the MSSM sector. Again in the 
case of AsyDM the relic density of Xi must be much smaller than the WMAP relic density for CDM. 
To this end we carry out an explicit analysis of the relic density within supergravity (SUGRA) grand 
unification |39j. As we show in Fig. |2]the relic density of Xi can be very small (see the next section 
for more detail), which allows the dominant component of the dark matter observed today to be the 
asymmetric dark matter. 



7 Detection of dark matter 

The AsyDM in the model we consider can interact with the standard model particles only via the 
Z' and p bosons which couple with the second and third generation leptons. Thus the scattering 
of AsyDM from nuclear targets will not produce any visible signals and the detection of AsyDM in 
direct detection experiments is difficult. However, it is interesting to investigate if the subdominant 
component of dark matter could still provide a detectable signature. We discuss this topic in fur- 
ther detail below. First we discuss the depletion of Xi dark matter to determine the regions of the 
parameter space where the relic density of Xi is a negligible fraction of the WMAP relic density for 
CDM, and is thus indeed a subdominant component of dark matter. Later we will investigate the 
possibility of detection of this subdominant component in direct detection experiments. Specifically 
we investigate two classes of models: the supergravity grand unified model (mSUGRA) with uni- 
versal boundary conditions on soft parameters at the GUT scale, and non-universal SUGRA model 



with non-universalities in the gaugino sector (see, e.g., 40 and the references there in 



For the mSUGRA case the parameter space investigated was: niQ < 10 TeV, mi/2 < 10 TeV, 
|^o/™o| < 10, 1 < tan/3 < 60, and > 0. For the case of SUGRA models with non-universalities 
in the gaugino sector the parameter space investigated was: Mi = m\j2 (1 -|- (5j) with the same 
ranges as in the mSUGRA case with < 1. After radiative breaking of the electroweak symmetry 
we collected roughly 31.4 million mSUGRA models and 25.6 million non-universal (NU) gaugino 
models. These models were then subjected to the experimental constraints which included the lim- 



its on sparticle masses from LEP 32 : m^^ > 81.9 GeV, m-± > 103.5 GeV, mr > 95.7 GeV, 
m^^ > 89 GeV, mg^ > 107 GeV, > 94 GeV, and mg > 308 GeV. as well as the recent bounds 
on the light CP even, SM-hke, Higgs from ATLAS and CMS, i.e 115 GeV < mf,o < 131 GeV (4ll[42 . 
More recently, ATLAS has constrained the SM-like Higgs to be between (117.5 — 118.5) GeV and 
(122.5 - 129) GeV and CMS has constrained the Higgs Mass to be between (115 - 127) GeV |43]. 
This constraint applied to the mSUGRA parameter space has recently been discussed in |44| . As 
discussed above if Xi < Xi* then the neutrahno would contribute to the relic density and for the 
AsyDM model to work we require that ^^ohl, to be less then 10% of the WMAP-7 limit [ij. 
Other constraints applied to the parameter points include the — 2 [33] constraint discussed in 



23 



0.008 



0.004 
0.002 



("1 1 v "\X71\ T A "P 
S , U.iXVViMAr^ 








'. - '*■ -..'■ - 


":>■ 


■A- 






•v",— 
f \ 


. '.' . 




v'. mSUGRA 







10 



10" 



10" 



10 ■ 



10" 



60 80 100 120 140 

Afy, (GeV) 



XENON- 100 



^ XENON-IT 

SupraCDMS-iT^ 

n = n^«hl/(nuMhl) 

mSUGRA 



40 60 



10 100 120 140 

Afj„ (GcV) 



126 
124 

3 

'o °>< 
120 ^ e 

1118 
116 



0.01 
0,008 
0.006 
0.004 
0.002 



O.lxWM-AP 



125 3 



\" NU gauf>;ino 




Q 
< 



60 80 100 120 140 

A/jo (GeV) 











10 




126 






10-" 




124 










122 






10-"'' 






















Nl20 




b 


10-" 






< 


X 














118 














10-=° 




116 
















in"== 



XENON-lOO 



* fei- XEN9N-rr__ 

SuperCDM'"' i 



7^ = nJofe2/(nDM/i 



40 60 80 100 120 140 

A/j,o (GcV) 




Figure 2: (color online) Top panel: An exhibition of the depletion of the MSSM neutralino dark 
matter below 10% of the WMAP relic density for cold dark matter. Parameter points are displayed 
by their light CP even Higgs mass and the yellow band corresponds to 10% of the WMAP-7 observed 
limit. Bottom panel: An exhibition of the neutralino-proton spin-independent cross section as a func- 
tion of the neutralino mass. To account for the reduced relic density of the neutralino component of 
dark matter the spin-independent cross section has been corrected by a factor TZ = JI^o/iq/ (J^dm^o) • 
The present experimental limits (solid line) f52] as well as the future projections (dashed lines) are 
shown (53{|54| . The left panel shows the parameter points of mSUGRA and the right panel shows the 
non-universal gaugino parameter points. All parameter points shown pass the general constraints. 



Section |5] and constraints from B-physics measurements 45 -47 which yield flavor constraints, i.e. 
(2.77 X 10~^) < Br {h ^ 57) < (4.37 x 10"'^) (where this branching ratio has the NNLO correc- 
tion j48|) and Br{Bs — t- /i^/x~) < 4.5 x 10~^. As done in 49 , we will refer to these constraints 



as the general constraints. These constraints were done by calculating the sparticle mass spectrum 
with SuSpect [50] and using MICROMEGAs [51] for the rehc density as well as for the indirect 
constraints. 

In Fig. [2j we exhibit the mSUGRA (left panels) and the NU gaugino (right panels) parameter 
points after applying the general constraints. In the top two panels we show the thermal relic 
density of the neutralino. As discussed previously, one finds that there is a significant region of the 
parameter space with a relic density much less than one tenth of the WMAP relic density. Thus the 
neutralino is indeed a subdominant component of dark matter. There are many more NU gaugino 
parameter points that satisfy the relic density compared to mSUGRA parameter points. This comes 
about because of coannihilation. Thus the non-universal case allows for the light chargino to lie close 
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Figure 3: Left: Leptonic final states in a collider where the — )• with / = fi, f^, r, Ur, 

final state arising from direct channel poles involving Z' . The Z' pole does not allow for a e~^e~ 
final state and thus the relative production cross section for fJ,^fi~ — >■ r^r~ vs — s- e^e~ can 
be used to detect the existence of a — Lr gauged boson. Right: A similar analysis is possible for 
p but its production is suppressed relative to Z' since it must be produced at the loop level. 



to the LSP, i.e. (m-± — m^o)/m^o <C 1, allowing for coannihilation to occur. The relevant question 

Xi Xi X-i 

then is if such a subdominant component can be detected in dark matter experiments. This is 
exhibited in the lower two panels of Fig. [2] where the corrected neutralino-proton spin-independent 
cross section, i.e TZ x o"|op where TZ = fi^o/ig/ (f^oM^o)' given as a function of the neutralino mass. 
For comparison we show the current XENON-100 bound |52| and the projected bounds in future 



experiments [53,54 . The important observation is that neutralino-proton spin-independent cross 
section is still detectable even when the neutralino is a subdominant component of dark matter with 
a relic density less than 10% of the WMAP relic density for CDM. 



8 Signatures at colliders 

The AsyDM models discussed above can produce a dramatic signature at a muon collider, see Fig. [3j 
which we now discussed (Signatures of a Z' boson in a gauged ?7(1)l^-Lt model at a muon collider 
were discussed in [22] but the analysis was only at the tree-level.). This signature arises from a Z' 
resonance. We note that Z' has no couplings with the first generation leptons and thus a process 
such as e^e~ — )• Z' — )• p,~^ ,t~^t~ is absent at the tree-level. This process can only arise at the 
loop level which, however, is suppressed relative to the tree. This explains why such a resonance has 
not been observed yet at an e~^e~ collider (see Appendix [c| . However, dramatic signals will arise at 
a muon collider where we will have processes of the type 

Since the final states contain no e~^e~ this would provide a smoking gun signature for the model. 
In Fig. g we exhibit the cross section a{fi'^fi — )• r+r ) for various values of gc when the AsyDM 
mass is taken to be 11.11 GeV and the Z' mass is 150 GeV. For comparison a[^^ jj,^ — )■ e^e~) 
is also plotted. One finds that the a{fi~^iJ,~ — >• r+r") exhibits a detectable Z' resonance and the 
cross section varies dramatically as a function of relative to a{fi~^fi~ — )• e~^e~) which is a rather 
smoothly falling function beyond the Z boson pole. In Appendix [C] it is shown that the loop 
contribution to /i^/i^ — s- e^e~ is suppressed and the Z' resonance is not discernible in this channel 
at a fJ,^fJ,~ collider. We note that there is a second overlapping resonance from a spin p state 
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Figure 4: (color online) An exhibition of the relative strength of the t~^t~ vs e'^e~ signal at a 
muon collider. The presence of a detectable Z' resonance in the ^'^ li~ — )■ t~^t~ channel provides a 
smoking gun signature for the gauged — L-t AsyDM model. A similar resonance is also present 
in the fi^^^ — s- fJ-^IJ'" channel while fJ.'^fi~ — )• e'^e~ cross section shows no such enhancement in the 
Z' region. 



where fj.'^fi — t- p — t- fl*fl — t- fi^fj, 1^ . However, the p resonance can only proceed at the loop level 
and is suppressed relative to the Z' pole. 



9 A gauged B — L model 

Next we discuss briefly the case where in the Stueckelberg extension we use U{1)b-l rather than 
U{1)l^-L^. Here we consider models with the right-handed neutrinos in order to gauge B — L and 
focus on the model K'2 with the B — L transfer interaction £asy = jvP^'^ i^^)'^ above the EWPT 
scale. By including three generations of right-handed neutrinos, the asymmetric dark matter mass 
is computed to be 6.06 GeV, c.f. Table [3j In this case, there are more experimental constraints to 
consider which include collider (i.e., LEP, Tevatron, LHC) constraints as well as precision constraints 
(i.e., the measurements of the p parameter, the T width). Specifically the LEP constraint gives 
Mz' /g'c ^ 6 TeV [55] for heavy gauge bosons. A stricter bound within a specific framework is given 
in ^56j where Mz' > 10 TeV. For lighter gauge bosons, as in the case of [9)|57], the UA2 cross section 
bound |58| is more stringent. Our analysis here is consistent with these constraints. Now, as in the 
— Lr case, the thermal symmetric contribution to the relic density from AsyDM must still be 
consistent with WMAP, i.e. it must be depleted to below 10% of the WMAP-7 value. An analysis of 
this is given in Fig. [sjfor the model Ag with 7 = (dashed line) and 7 = 70 = 1-3 x lO"^*^ (solid line). 
Here one finds that the symmetric component of AsyDM can satisfy our WMAP-7 constraint for a 
range of Z' masses. If one wishes to keep the Z' mass in the T region, i.e. ~ 10 GeV, then a fine-tuned 
value of g'fj is needed as seen in Fig. [sjto be consistent with the constraints on Mz'/g^- Heavier Z' 
masses would have difficulty satisfying both the collider and WMAP-7 constraints discussed above. 
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Figure 5: (color online) A display of the thermal relic density of ^0 as a function of Mz' for the model 
with a gauged B — L for different couplings with 7 = (dashed line) and 7 = 7o = 1-3 x 10"^*^ (solid 
line). It is seen that resonant annihilation of thermal dark matter via the Z' pole allows the relic 
density of this component to be reduced to below 10% of the WMAP value for values of Z' around 
twice the mass of the dark particle. 



10 Conclusion 

In this work we have proposed models of asymmetric dark matter in the framework of Stueckel- 
berg extensions of the SM, the 2HD and the MSSM. Several candidate models for asymmetric dark 
matter were discussed using a variety of operators constructed out of standard model fields which 
carry a non- vanishing B — L quantum number which is transferred to the dark matter sector at high 
temperatures in the early universe consistent with sphaleron interactions which preserve B — L. The 
analysis was done both for models where the interaction temperature at which the B — L transfer 
takes place above the electroweak phase transition scale as well as below this scale. The details of 
the B — L transfer determine the mass of the asymmetric dark particle. A master formula was given 
which generates the asymmetric dark matter mass for a variety of models discussed in the text and 
allows one to discuss new possibilities. Specific models are discussed, including those anchored in 
the standard model, the two Higgs doublet model as well as the minimal supersymmetric standard 
model, with or without the right-handed neutrinos. 

A central ingredient in a successful asymmetric dark matter model, and an explanation of cosmic 
coincidence, is an exhibition of an efficient mechanism for the annihilation of the symmetric dark 
matter component which is produced by thermal processes. We accomplish this using a Stueck- 
elberg extension of the standard model and of the MSSM. The Stueckelberg extension of the SM 
is particularly simple and appealing in which, aside from the dark matter field, there is just one 
more field, a gauge boson (Z') of a — Lr gauge symmetry, which gains mass via the Stueck- 
elberg mechanism. The symmetric dark matter produced by thermal processes is depleted via 
resonant annihilation from the exchange of a Z' using a Breit-Wigner pole. This is perhaps the 
simplest asymmetric dark matter model, in that there are no extra Higgs fields that appear in the 
model. Moreover, the extra U{1) gauge symmetry in our models forbids the dangerous Majorana 
mass terms which would generate oscillations of the dark particles and their anti-particles which 
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could washout the asymmetry. 

A supersymmetric extension of this model is also given where it is shown that in addition to the 
Z' boson, there is also a spin boson field along with two additional neutralino states arising from 
the Stueckelberg gaugino sector. It is shown that the p has a rapid decay and does not participate 
in the dark matter analysis. In the MSSM extension, there is an extra complication, in that, with 
R parity one can have a stable neutralino which is a possible dark matter candidate, and it must be 
shown that it is also depleted so does not compete with the asymmetric dark matter candidate. In 
the analysis presented in this work it was shown that there exists a significant part of the parameter 
space of mSUGRA where the relic density arising from the neutralino was less than one tenth of the 
WMAP relic density and thus the neutralino is a subdominant component of dark matter. Inter- 
estingly, however, it was shown that the subdominant neutralino is still accessible at future direct 
detection experiments such as SuperCDMS and XENON-100. It was shown that definitive tests of 
the model can come from a muon collider where one can produce the Z' which decays only into /i's 
and r's (and /x, r neutrinos). 

We also discussed a gauged B — L Stueckelberg model. Again in this model the symmetric 
dark matter can be efficiently depleted by annihilation near a Breit-Wigner Z' pole. Thus within 
the Stueckelberg extensions there exist several possibilities for explaining cosmic coincidence. The 
dominant dark matter in all these models will be light and lies in the range of 1 — 16 GeV. 
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Appendix 

A Master formula for computing the asymmetric dark matter mass 

We have discussed various Models A-F and subcases such as Ai-Ag etc, and also models with 
right-handed Dirac neutrinos, which will be discussed at the end of this section. We discuss now 
a master formula which allows one to take some particles in or out of thermal equilibrium. Such 
a formula would generally be useful before SU{2)l breaking, i.e., T > Tewpt, when some of the 
super-particles are suppressed in the plasma while others are not. This would allow us to discuss the 
Models A,D,E,F in a unified way and also allow us to generate new models where some other sets 
of super-particles are taken out of the relativistic plasma in the early Universe. However, it is not 
useful to discuss such a formula below the electroweak phase transition scale since the current experi- 
mental data indicates the sparticles to be heavy and not below the electroweak phase transition scale. 



In obtaining the master formula, we assume: (1) In supersymmetric cases, all particles in a 
supermultiplet have the same chemical potential; (2) A given particle type in different generations 
has the same chemical potential, e.g., fi^ = fJ-s = fJ'b] (3) All the additional Higgs doublets have 
the same chemical potential as the standard model Higgs p,H- Following the discussion of Sec- 
tion 3.1 for all the fields in the plasma, we have the chemical potential constraints as before, i.e. 
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fiH = Hl — fJ-e = fJ-q — fJ-d = f^u — fJ'q from Yukawa couplings, 3fiq + fj,L = fr'om sphaleron processes, 
and Y = from the neutrality condition. 

(i) 

It is useful to introduce the temperature-dependent coefficients Ca for the matter fields in the 
plasma. We define Ca^ = Ca\f) + Ca\b), where Ca\f ) counts the contribution of i**^ generation 
particle a (with mass m^) which is fermionic and Ca\b) counts the contribution of its super-partner 
a (with mass rria) which is bosonic, where Ca{f) and Ca{b) are given by Eq. Q. We note that in 
the limit when rua can be neglected, one has a weakly interacting plasma so that Ca{f) = 1 and 
Ca{b) = 2. Thus we have, for T ^ rria, c~a = 1 -|- 2 = 3; for ^ T ^ ttIq,, c~a =1 + = 1; and 
for T <^ nia, Ca^ = 0. For the Higgs doublets, we have ch = cnib) = 2 in the non-supersymmetric 
case, and ch = cnib) + cnif) = 3 for the supersymmetric case, and Xh counts the number of Higgs 
doublets in the model. 

We can then rewrite the hypercharge neutrality condition as 

2(c« + + cf + 4(c(i) + c(2) + c(3))^. - 2(c« + + c?V. 

-2(4') + + - 2(c(') + c(2) + cf));,, + 2cHXHf^H = . (114) 

Further defining Ca = '^c? = Ca ) + c'a^ + ci\ and together with Eqs. (p0l)-(p2|, we obtain 



and 



B-L 



= 6C. + 3C, + 3a + 3c^AH "^^-^^^ ' ^'''^ 



: - C„(3C, + + 9Cl + 2Ce) + Cd(3C7g + 9Cl + 8Ce) + Ce{Cq + 3Cl) 
+ Ci^ Ah {2Cg + Cu + Cd + 6Cl + 3Ce)]fiL / (6C„ + 3Cd + 3Ce + 3chAh) 
: -AA/iL/(6C„ + 3C7rf + 3C7e + 3c/^A/^) . (116) 



From Eqs. ([8|) and (28), we find the master formula for computing the dark matter mass 

"^DM ^ ^ + 3Cz. + 6Ce + 3c^,A^^ ' -gg_^ ' 97 ^""^ ' ^^^^^ 



where k is the parameter indicating the dark matter type: for the non-supersymmetric case, k = 1 
for fermionic dark matter, and k = 2 for bosonic dark matter; for the supersymmetric case, k = 3. 
Q%_l is the [B — L)-charge of the operator Oasy in the B — L transfer interaction Eq. ([2|. The 
operators that carry Q^_i = — 1 are: Oasy = LH, LLe'^, LqdF, u'^cFdP. Higher dimensional operators 
with a larger value of {B — L)-charge can be constructed from them. 

We can extract the results for Models A,D,E,F from this master formula: 

1. Model A: For matter fields we have c'a = Ca\f) = 1 (so that Ca = 3); ch = cnib) = 2 and 



Xff = 1. For fermionic dark matter we take k = 1 and we recover Eq. (35). 



2. Model D: Here all the parameters are the same as in Model A except that Xh = 2. Setting 



K = 1 we recover Eq. ( 37 ) . 
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3. Model E: For matter fields, S = S{f) + c^a{b) = 3 (so that Ca = 9); ch = CH{b)+CH{f) = 3 
and Xh = 2. Since this is a supersymmetric case, k = 3 and we recover Eq. ( |69[ ). 

4. Model F: For matter fields, ^ = Ca^ = 1, Ca^ = 3 (so that Ca = 5); ch = 3, \h = 2. Taking 



K = 3 we recover Eq. ( 69 ) 



For Models A',D',E',F' which include three generations of right-handed Dirac neutrinos in the 
thermal bath, the master formula reads, 

C, + 8C. + 2Q + 3C, + 6a + 3cHAH + J " ^QoT, ' ^ ' ^'''^ 
where Ci,j^ = ci^ + c^jl + c^u^. By taking (7,^^ = 3,3,9,5, we can recover the dark matter mass 



formulas for Models A',D',E',F' (Eqs. p9} and (73)). 



We can also obtain the dark matter mass of other models from the master formula Eq. (117) or 



Eq. (118) by varying the temperature where certain heavy particles are Boltzmann suppressed in 



the thermal bath. 



B Decay of the p 



Here we compute the decay of the p. From Eq. (112) one finds that p couples to smuons, staus, 
muon sneutrino, and tau sneutrino. This means that the p decay has i^^z^^, t~^t~, VtI^t final 

states which arise via the exchange of neutralinos and charginos in the loops (a generic diagram is 
shown in Fig. [g]). The amplitude of the generic diagram reads, 

= -^^P^^-^^^^^. / (2^-(^ )(fc2_^2)(^,2_^2)((^_^)2_^|j-(^') ' (119) 

where k' = q — k, mi,mj are the masses of the sleptons, and m^^ is the mass of the neutralino or of 
the chargino in the loop, while gpij, C^i are the couplings. Our purpose here is to estimate the size 
of the lifetime and to that end it is sufficient to estimate the contribution for one set of diagrams. 
Thus we consider the decay of the p to final states p^ p~ via the exchange of neutralinos. In this 
case we will have the exchange of smuons and neutralinos in the loop. Further, we will ignore the 
mixing between the left and the right chiral smuons so that the mixing angle Or in Eq. (112) can 
be set to zero. In this circumstance the off-diagonal term involving two smuons in the loop 
contribute and the relevant loop integral takes the form 



does not 



1 



2 

(12 2\(Mi 2\t(i ^2 2\- dxdydz6{x + y + z - 1)—^ 
where Dik = — Aj^ -|- ie in which I = k — {yq + zp) and 

,2 _„._2 , ^2 , / 2 „\_2 



(120) 



Aik = (1 - z)mi - xynip + zm^^ + {z - z)mp . (121) 

The masses in the loops are much larger than the muon mass and thus the muon mass can be 
ignored. The integration on I gives 

^M = f axdydz6ix + y + z- l)^M^x^ . (122) 

(47r)^ 7o ^ik 
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k' 



u{p') 



k-p 



v{p) 



Figure 6: A generic diagram showing the decay of the p to one of the final states which could be 
fj,^p^ , VfjiT^fji, T~^T~ , VrVr via exchange of sleptons, charginos and neutralinos at one loop. 



Further, an approximate evaluation of integration on the Feynman parameters gives 

iM 



-}M9^u{p')'^v{p) , 



(4^)2 



(123) 



under the assumption rn^^/mf <^ 1 and m^/m? <^ 1. The decay width of p — )■ fx^fi is then given 
by 



dr 



2m. 



d'^p 



{2tt)'^2E + J (27r)32K 



\EiM\ 
Sirnin 



(124) 



Next we note that gpu = —gp22 = dcQcm-p and thus 



2 {gcQc)^m^ 



167r4 



6 2 



k=l 1=1 



2 



(125) 



A numerical estimate using Eqs. (124) and (125) and the inputs mi = 1 TeV, m2 ^ mi, mp = 
100 GeV, the lightest neutralino mass of 50 GeV gives Tp = h/F 10~^^^^ s. Thus the decay of the 
p is very rapid. 



C Z' exchange contribution to fi~^fi e^e at loop level 

At a muon collider, e~^e~ final states can be created via photon exchange and via a Z boson exchange. 
Since the Z' has no direct coupling with the first generation leptons, there is no tree-level Z' exchange 
contribution to e~^e~ final states. However, at the loop level a Z' exchange can make a contribution 
where the second and third generation leptons are exchanged in the loop as shown in Fig. [7] We 
now compute this contribution to determine its size. Thus we consider a p-^p-^ — t- e~^e~ process 
with a Z' exchange via the second and third generation leptons loops as shown in Fig. [7] In this 
case the contribution to the scattering amplitude is 

= iMyZ' +iMzz' 

= ^(pO(^gcQ.7n^(p) "'^^^"'"t^"^'^ ^n7z0^n(fc)(-.e7-).;(fcO 
+ v{p ){i:9cQpT)u{p) ^ 2 V^zz') 
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Figure 7: Z' — 7 and Z' — Z exchange via [i^ [i , i^^f^, t~^t , l'tI't loops. 



2 , u{k)^i^igv+gAl'Hk'), (126) 
— cos t)w 

where is the U{1)L^-Lr charge for muon, gv = ^(r3)i + sin^ 9w, 9A = —5(^3)1,, and the vacuum 
polarization tensors iH^^, and iH^'^, are the sums of the contributions from /z"*"//", z^^z^/^, t~^t~ , Vt-Vt 
loops 

iii^p = in^ + m;;^'; + iw/ + zn;;^ . (127) 

First we focus on ^Hj^^, ^ which is the muon loop contribution to the Z' — 7 exchange. It reads 



Kz'u = -ii^9cQ,){-^e) I ^tr 



^i{f + m^) i{f' + m^) 



7Z>- ^2^^'^^^^ (2^)4^1^ (r2_^2)T (^,2_^2) 



(47r) 



{q^g^P-q^qn [ dx x(l - (- - logA^ - 7 + log(47r) + 0(e)) , (128) 
JO e 

where = — x(l — x)(;^, and in the last step we use the dimensional regularization. The 
expression of iS^^' r 'ii^^i's from iS^^, ^ by only the Qr factor, and it takes the form 

^Kz',r = ^^^^i'l'^a"' - 'f^') [ dx x(l - x) - logA^ - 7 + log(4^) + 0{e)) . (129) 

Summing over these two terms, we find a dramatic cancellation of the divergence in the loop due to 
Qp = —Qt = 1) making the loop finite so that 



where 

dxx(l — xjlog- — = / dxxfl — xjlog — 7^ (131) 

A^ Jo - - x)q 

One can also obtain the neutrino exchange contributions from the above by setting the fermion 
masses to zero in the equation above (assuming neutrinos to be massless) which gives a vanishing 
contribution. 



32 



Now we want to compare the contribution of the Z' — 7 exchange loop diagram with the tree-level 
process fJ,'^fJ,~ ^ ^ e^e~ , whose amplitude reads 



iM^ = v(p')(^-iejt')u{p)^§^u{k)i-ieY)v{k') . (132) 
With some manipulation we find 



Hi 



.2 



^•^72' = - (gj2 • g2 _ ^2 ^ X iA^^ = / X iM^ . (133) 
Thus, the total squared amplitudes involving a photon can be written as 

\iM^ + iMjZ'l^ = |1 + /P X |«A^7p 

= {l + f + r + fnx\iM^\'- (134) 

Our numerical analysis shows that (/ + /*+ //*) is smaller than ~ 10~^ and thus the loop makes 
only a tiny contribution to the total cross section in this case. The analysis of Z' — Z exchange is 
similar and gives a very small value. Thus we conclude that a Z' peak will not be visible in the 
e+e~ process at a muon collider. The above analysis also exhibits why a Z' in this model 
would not be visible in an e+e~ machine. 
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